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ABSTRACT 

We study the dynamics of the 3-D three-body problem of a small body moving under the 
attractions of a star and a giant planet which orbits the star on a much wider and elliptic or- 
bit. In particular, we focus on the influence of an eccentric orbit of the outer perturber on the 
dynamics of a small highly inclined inner body. Our analytical study of the secular pertur- 
bations relies on the classical octupole hamiltonian expansion (third-order theory in the ratio 
of the semi-major axes), as third-order terms are needed to consider the secular variations of 
the outer perturber and potential secular resonances between the arguments of the pericenter 
and/or longitudes of the node of both bodies. Short-period averaging and node reduction (by 
adoption of the Laplace plane reference frame) reduce the problem to two degrees of freedom. 
The four-dimensional dynamics is analyzed through representative planes which identify the 
main equilibria of the problem. As in the circular problem (i.e. perturber on a circular or- 
bit), the "Kozai-bifurcated" equilibria play a major role in the dynamics of an inner body on 
quasi-circular orbit: its eccentricity variations are very limited for mutual inclination between 
the orbital planes smaller than ~ 40°, while they become large and chaotic for higher mutual 
inclination. Particular at tention is also gi ven to a region around 35° of mutual inclination, de- 
tected numerically by .Funk et al.l (1201 ll) and consisting of long-time stable and particularly 
low eccentric orbits of the small body. Using a 12th-order Hamiltonian expansion in eccen- 
tricities and inclinations, in particular its action-angle formulation obtained by Lie transforms 
inlEibert & Henrard (2008), we show that this region presents an equality of two fundamental 
frequencies and can be regarded as a secular resonance. Our results also apply to binary star 
systems where a planet is revolving around one of the two stars. 

Key words: planetary systems, celestial mechanics, methods: analytical, planets and satel- 
lites: dynamical evolution and stability, binaries: general 



1 INTRODUCTION 

The (inner) Lidov-Kozai mechanism (Kozai 1962; Lidov 1962) is a 
well-known secular resonance of the restricted three-body problem, 
which can be reduced two degrees of freedom after short-period 
averaging and node reduction (see for instance Malige et al. 2002). 
Kozai (1962) showed that a highly inclined asteroid perturbated 
by Jupiter periodically exchanges its eccentricity and inclination. 
Its analytical theory relied on the assumption that Jupiter's orbit 
is circular, so that the problem is integrable. Since its discovery, 
the Kozai resonance has found numerous applications in studies of 
planetary and stellar systems. 

Recently, analytical studies (e.g. Michtchenko et al. 2006, Lib- 
ert & Henrard 2007) have shown the possibility that extrasolar plan- 
etary systems can be in a long-term stable highly non-coplanar con- 
figuration, sometimes due to a secular Kozai-type phase-protection 
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mechanism. For instance, Libert & Tsiganis (2009) found that 
V Andromedae, HD 12661, HD 169830 and HD 74156 extraso- 
lar two-planet systems have orbital parameters compatible with a 
Kozai-resonant state, if their (unknown) mutual inclination is at 
least 45°. 

The Kozai dynamical phenomenon is also well-known in stud- 
ies of binary systems (e.g. Innanen et al. 1997, Wu & Murray 2003, 
Fabrycky & Tremaine 2007), in particular in the S-type configura- 
tion (a planet revolves around one of the primaries) where the orbit 
of a highly inclined planet can undergo large amplitude oscilla- 
tions of its eccentricity. A similar Kozai-resonant evolution can be 
observed in the C-type configu ration, where the planet is in orbit 
around the binary (e.g. Migasze wski & GozdziewskiluOllh . 

Concerning the planetary three-body problem, the discovery 
of giant extrasolar planets on eccentric orbit raises the question of 
the influence of their eccentricity on potential asteroids and Earth- 
mass companions on inclined orbit. In a preliminary numerical 
study of Funk et al. (201 1), the long-term stability of inclined fic- 
titious Earth-mass planets in the habitable zone of extrasolar giant 
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planets discovered so far is analyzed. They have realized a paramet- 
ric analysis, where several giant planet's eccentricities are consid- 
ered, while the Earth-like body is initially on a circular orbit, closer 
to the star than the gas giant. Their simulations show that, for dis- 
tant orbits, test-planets below a critical inclination of approximately 
40° are in a stable configuration with gas giants on either circular 
(i.e. well-known result associated with the Kozai mechanism) or el- 
liptic orbit. Furthermore, for gas giant on eccentric orbit, the small 
companion exhibits non-negligeable variations in eccentricity, ex- 
cept for a region around 35° of mutual inclination of the orbital 
planes, consisting of long-time stable and low eccentric orbits of 
the Earth-like body. So the influence of the eccentricity of the per- 
turber on the dynamics of inclined Earth-mass planets seems to be 
significant and deserve to be studied in more detail with dynamical 
tools. This is the goal of the present contribution. 

In this work, we focus on the 3-D three-body problem of a 
small body moving under the attractions of a star and a giant planet 
which orbits the star on a much wider and elliptic orbit. Our an- 
alytical study of the secular perturbations relies on the classical 
octupole hamiltonian expansion (third-order theory in the ratio of 
the semi-major axes), widely used in planetary and stellar sys- 
tems (e.g. Ford et al. 2000, Lee & Peale 2003, Migaszewski & 
Gozdziewski 2011). Actually, third-order terms are needed to in- 
troduce the secular variations of the eccentricity of the perturber. 
Indeed, the second-order quadrupole approximation does not de- 
pend on the argument of the pericenter of the perturber, whose 
eccentricity is thus an integral of motion (e .g. lHarri ngton*4969. 
Lidoy&Zigiir]! 19761 iFerrer & Osacaijll994l and|Farago & Laskar 
20101) . The third-order terms introduce qualitative changes in the 
dynamics and can explain the aformentioned dynamical features 
observed by Funk et al. (2011), as we will show in this work. Let 
us note that, even if the octupole development is an analytical ex- 
pansion of the three-body problem whatever their masses, we only 
focus on planetary systems with a small value of the inner body's 
mass hereafter. This problem is sometimes called the reduced spa- 
tial three-body problem. Since we consider the (very limited) ef- 
fect of the small mass on its companion, the secular evolution of 
the outer body is considered, and so are the potential secular reso- 
nances between the arguments of the pericenter and/or longitudes 
of the node of both bodies. 

The paper is organized as follows. In Section 2, the octupole 
Hamiltonian formulation is recalled. Section 3 analyzes the four- 
dimensional secular dynamics of the elliptic spatial three-body 
problem, by means of 2-D geometric representations called rep- 
resentative planes. Section 4 focusses on the dynamical feature 
around 35° of mutual inclination of the orbital planes, described 
in Funk et al. (201 1). Finally our results are summarized in Section 
5. 



2 OCTUPOLE HAMILTONL^N FORMULATION 

Let us consider a system consisting of an inner small body (mi) 
and an outer giant planet (m 2) orbiting a star (mp) ( also called in- 
ner three-body problem, see F arago & Laskarj 120101) . Due to their 
masses, the inner body will be named the perturbed body and the 
outer one the perturber in the following. We focus on the spatial (or 
3-D) problem where both planetary orbits are mutually inclined. 
In the Solar System, this configuration corresponds for instance 
to an asteroid perturbated by Jupiter Let us note that studies on 
the secular evolution of asteroids are mostl y realized un der the as- 
sumption that Jupiter's orbit is circular (e.g. lKozaill962h . Since the 



discovery of extrasolar systems, inclined test particles, represent- 
ing Earth-mass planets with weak gravitational effects on a system 
composed of a star and a gas giant, are another application of the 
spatial problem. However, as many giant extrasolar planets have 
eccentric orbit, one may wonder the influence of the eccentricity of 
such a Jupiter-like planet on its Earth-mass companion(s). To ad- 
dress this question, we consider in the following that perturber is 
on an eccentric orbit. 

The spatial model of the three-body problem can be described 
using the canonical hel iocentric formulation (see .Poincare.. 1 89q , 
iLaskar & Robutelll995t) : 



^ I ?m'. 



J=l 



pj G(rao + mj)m'. 



m\m2 

l|ri -r2ll 



+ ^^. (1) 

mo 



where r, are the position vectors of m, relative to the star, p, are their 
conjugate momenta relative to the barycenter of the three-body sys- 
tem, and m'- = (1/mo-l- 1/m;)"' are the reduced masses. Let us recall 
that the heliocentric velocities r, will not follow the direction given 
by p,, and thus the ellipses are not tangent to the real trajectory. 
The first term of the expansion is the sum of the keplerian motions 
of the two planets. The perturbation of this integrable part, repre- 
senting the mutual interactions between the planets, consists of the 
direct part and the indirect part respectively. 

A set of canonical variables is formed by use of the classical 
Delaunay's elements: 



Sj 



■ Mj, Lj = m'j -^G(mo + mj)aj, 

■ (Dj, Gj = Lj -y'l - e/. 



(2) 



hj = ilj, Hj = Gj cos ij. 



where aj denote semi-major axes of the planets, Cj eccentricities, 
ij inclinations, ojj arguments of the pericenter, Qj longitudes of as- 
cending nodes, and Mj mean anomalies, all being canonical helio- 
centric elements. 

As we are interested in the long-term dynamics and assuming 
that we are not close to a mean motion resonance, we can aver- 
age (to first order in the mass ratio) the Hamiltonian function over 
the fast angles, namely the mean anomalies M, (Deprit 1969). It 
means that the averaged Hamiltonian "K does not depend on the 
mean anomalies, then the conjugate momenta L, are constants in 
the secular problem and so are the semi-major axes. So it results in 
a four degree of freedom formulation of the Hamiltonian function. 

To average the indirect part of the disturbing function, we 
compute ji^ Jo" Jo" Pi ■ Pj dMjdMj = SijO^n^, where p„ Pj are 
canonical heliocentric velocities related to the canonical heliocen- 
tric elements. 

Concerning the direct part, we use the traditional expansion in 
Legendre polynomials, assuming that ri « r2: 



1 °° / \" 
Hop = -Gmim2—}i — \ P„(cosS), 



(3) 



where 5 is the angle between vectors ri and r^. We choose to per- 
form the development for all P„ with n < 3. This well-known devel- 
opment, limited to order 3 in th e semi-major axes ratio a = fli/fl2, is 
c ailed octupole theory (see e.g.lFord et al.ll2000l iLee & Pealell2003l 
and Migaszewski & Gozdzi ewskiluOI ih . For the sake of complete- 
ness, we present hereafter the technical details of the calculations. 
Practically, the term cos 5 = ri • r2/(rir2) can be expressed in 
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the following way - see lPradd ( I2005I) for a similar 2-D calculation: 

(4) 
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where i?i (i'l , Wi , fii ) and /?2('25 ^2, n2) are Eulerian rotations of the 
orbital reference frames of the mass mi and mi respectively, and a 
and j6 have rather simple expressions: 

Q- = i?l.l cos/2 +i?i,2 sin/2, y6 = i?2,l cos/2 + /?2,2 sin/2, (7) 

Rjj being the element in the ith row and 7th column of the matrix R. 
The averaging of the Hamiltonian lO over the short-period terms, 



< Hop >Mi,M2 ■ 



1 r^2n r>2jl 

(2^ Jo Jo ^' 



opdMidM2 
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is realized, in practice, over the eccentric anomaly Ei of the per- 
turbed body and over the true anomaly /2 of the outer perturbing 
body. The intermediate results after the first averaging are the fol- 
lowing: 
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+ 15(<5^+y3^)- 12], (12) 



where Hopt means the term of the direct part (|3) associated to the 
ith Legendre polynomial. For the averaging over the true anomaly 
of the outer body, we first replace 5- and fi by their values (see Eq. 
Q), and obtain the following first terms of the secondly averaged 
Hamiltonian: 



< Hnpn > 



Gmi mi 



DPO ''Ml. Ml 



ai 



< Hnpi > 



DPI ^Mi.M, 



0, 



(13) 
(14) 



the terms < Hopi >mi,Mi ^nd < H^pi >miMi being too long to 
be exposed here. As we can see, the first two terms are constant in 
the secular problem and do not contribute to the averaged Hamil- 
tonian TC. Let us note that an alternative development of the secu- 
lar expansion using Hansen coefficients can be found in Laskar & 
Boue (2010). 

To simplify the formulation of the averaged Hamiltonian, the 
Jacobi's redu ction, also known as the elimination of the nodes 
l |jacobilll842l) . allows us to reduce the expansion to a two degree 
of freedom function only. The reduction is based on the existence 
of additional integrals of motion, namely the invariance of the total 
angular momentum, C, both in norm and in direction. The constant 
direction of the vector C defines an invariant plane perpendicular to 
it. This plane is known as the invariant Laplace plane. The choice 
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of this plane as reference plane implies the following relations: 

^1-^2 = ±180° (15) 



Gi cos /] + G2 cos 12 = C 

Gi sin ii - G2 sin !2 = 



(16) 
(17) 



with C being the norm of the total angular momentum. Let us note 
that, in our study, the invariant Laplace plane coincides almost with 
the Jupiter-like planet's orbital plane, since its inclination relative 
to the invariant plane, i2, is of the order of m\ jm'^ by relation Ml\ . 
Another quantity, related to the total angular momentum, is 
frequently used. This is known as the angular momentum deficit 
(Laskar 1997): 



AMD = ^ Lj{l - J\ - e] cos ij) 



■ Li+ U-C. 



(18) 



Finally, we present the octupole expansion of the Hamilto- 
nian O averaged over the short-period terms and expressed in the 
invariant Laplace plane, using the succinct formulation introduced 
bv lFord et al.l ( l2000h : 

^- --^ [(2-3(l-(|f))(3cos^.-l) 

+ 15(l-(fi)')(l-cos2/)cos2gi] 

+ " ^^^^ Vi-lir) Vi-(ir) 

[A(-cosgi cosg2 - cosisingi sing2) 

+ 10(2^1 cos ('(l - cos^ /) singi sing2 

(19) 
where 



2G1G2 



S = 7-5(fif-7(l-(fi)')cos2g, 



^7-3(irf 



1(1 -COSMOS, 



i = i\ + i2 being the mutual inclination. The equations of motion 
associated to Hamiltonian <19l l are 



cfK 
dGi' 



G, = - 






(20) 



One has to keep in mind that such an approach is limited to 
small values of the semi-major axes ratio, namely o- < 0. 1. To con- 
sider larger values of the ratio, a development to higher order is 
needed, as done bv lKozail (1962). 

For a Jupiter-like planet on a circular orbit (G2 = L2), the for- 
mulation ( |19t simplifies to the quadrupole approximation (second- 
order terms in a). Then the secular Hamiltonian does not depend 
on g2, and the norm of the associated momentum G2 is an in- 
tegral of motion, which means that the eccentricity of the outer 
body is constant in this formulation. As a result, the problem is 
integra ble and th is approximation is studied in many papers ( e.g. 
Harrin gton 19691 [Lidoy & Ziglin|[l97^ iFerrer & 0sacarj|l994l and 
Farago & Laskar 2010, or in the artificial satellite context e.g. 
Lidov 1962, Russell & Brinckerhoff 2009 and Delsate et al. 201(i. 

In the present work, no assumption on the eccentricity of the 
Jupiter-like planet is considered. The variation of the eccentricity of 
the perturber is introduced through the octupole terms (third-order 
terms in a). As these terms depend on the variable g2, the prob- 
lem can not be reduced to one degree of freedom, and it induces 
qualitative changes on its dynamics, as shown in the next section. 
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3 GEOMETRIC REPRESENTATION OF THE 
DYNAMICS 

In this section, we study tlie dynamics of the two degree of free- 
dom Hamiltonian l |19l l by means of 2 -D geometric representa- 
tions, called representa tive planes (see iMichtchenko et alJIZOOq . 
iLibert & Heiu'ardll2007h . The idea is to choose a 2-D plane of ini- 
tial conditions which is suitable for the analysis of the stationary 
solutions of the secular two degree of freedom problem. This plane 
should be representative in the sense that we aim to find a plane 
such that it contains the initial conditions of orbits representative of 
each class of orbits. 

Such a plane can be obtained by fixing g\ and gi to values that 
verify the conditions 



5^ 



and Gi 



ffK_ 
dgi 



0, 



(21) 



i.e. according to the symmetries of the secular 3-D problem, 
(2gi,Aai = gi - g2) = (0,0),(0,7r),(7r,0) and (;r,;r). Indeed, the 
secular Hamiltonian function can be developed in Fourier series of 
cosinus terms, whose generic argument is 



■ jigi + hgi + kAii, 



(22) 



where ji and J2 are of the same parity (ji,J2,k are integers), and 
Ail = n after node reduction. As a result, conditions (I21> are veri- 
fied when sin0 = 0, i.e. (2g,,Aw) = (0, 0), (0, ;r), (;r, 0) and (n,n). 
These four pairs of angles define four distinct quarters of the repre- 
sentative plane. 

In the fo llowing, we choose the geometric r epresentation in- 
troduced by .Migaszewski & Gozdziewskil ( 1201 ih , and defined as 
X = eicoslgi and )» = £2 cos Aw with sin 2gi = sin Aw = 0. On 
this representative plane, the level curves of constant Hamiltonian 
are plotted for given values of AMD, a = ail 02 and /j = m[/in2. 
The boundaiy of permitted motion is defined as i = 0, 180° here- 
after. Let us recall that the eccentricities and mutual inclination are 
related through the integral of AMD. 

We insist on the fact that the representative plane is not a 
phase space or a surface of section. However, all orbits have to 
cross the representative plane (i.e. pass through the conditions 
sin2gi = sin Aw = whatever the behavior of the angles 2gi 
and Aw), and the points of intersection have to follow a con- 
stant energy curve. As the extremal va lues of the e ccentricities 
are reached when sin 2gi = sin Aw = 1 Michtchenko et al. 2006 . 



iLibert & Henrardll2008i iLibert & Tsiganisll2009l) , a quasi-periodic 
solution intersects the representative plane at four points on the 
same energy level. A stationary solution, fulfilling the two addi- 
tional conditions g, = &K/dGj = (i = 1,2), appears as a fixed 
point on the plane, while a periodic solution for which an angle is 
fixed has only two points of intersection. Orbits of chaotic motion 
intersect it at an arbitrary number of points. Depending on the loca- 
tion of these intersection points on the four quadrants of the plane, 
the behavior of the angles can also be deduced, as well as an esti- 
mation of the ranges of eccentricity variations, as will be shown in 
the following examples. 

The mass ratio fi is fixed to 10""* in the following. Concern- 
ing the semi-major axis ratio, the use of the octupole terms limits 
the width of semi-major axes ratios that can be considered; to en- 
sure the validity of our appr o ach, w e choose a = 0.05. Indeed, 
iMigaszewski & Gozdziewskil ( 120 llh have shown that, for hierar- 
chical systems, the octupole formulation is very precise and higher- 
order contributions do not distort the structure of the Hamiltonian 
curves of the representative plane. 



Two initial configurations of the three-planet system are ex- 
amined in this work. In Section 13.11 the outer giant perturber is 
considered on a nearly circular orbit, while the influence of a highly 
elliptic orbit of the perturber is analyzed in Section l3!2l 

3.1 Nearly circular orbit of tlie perturber 

When the perturbing body is on a circular orbit, the octupole 
formulation coincides with the quadrupole approximation. Since 
the quadrupole approach '7C,„„rf(ei, wi) is integrable, its dynamics 
can be represented on the phase space (ei cos Wi,ei sinwi) (see 
e.g. iThomas & Morbidellil 1 19961 for more detail). However, for 
nearly circular orbit of the perturber, the expansion (1191 is four- 
dimensional and a first picture of the dynamics consists in the plot 
of the level curves of constant Hamiltonian in the aforementioned 
representative plane. This representation is given in Figure [T] for 
several values of AMD: 8.09 x 10"', 4.79 x 10"^ 1.11 x 10"', 
1.80 X 10"^ 2.683 X 10"' and 2.7 x 10"'. For all these values ex- 
pect the last one, the maximal mutual inclination between the two 
orbital planes, z„,„^, is reached at the origin (ei = 62 = 0), while the 
border enclosing the possible dynamics of the problem represents 
the coplanar case (; = 0). The five first AMD values considered 
here coiTespond to /„„„ of 20°, 50°, 80°, 110° and 180° respec- 
tively. In the bottom right panel of Figure[T] the region of permitted 
motion separates into two parts which are bordered by the curves 
i = 0° (higher absolute values of 62) and ; = 180° (smaller absolute 
values of £2)- 

As explained hereabove, the structure of the geometric repre- 
sentation reveals the equilibria of the problem. For small mutual 
inclination (!„,„j. = 20°, see Figure[T]top left), circular orbit of the 
inner body corresponds to a stable equilibrium and no variation 
in eccentricity is possible. For larger inclinations (/„,„ = 50° and 
80°, top right and middle left panels respectively), the point e^ = 
becomes an unstable equilibrium, and a separatrix divides the left 
panel of the representative plane (where 2gi = n): the closed re- 
gion is characterized by the libration of g, around 90° or 270° and 
the region outside the separatrix by the circulation of this angle. 
The two stable equilibria (at gi = 90° and gi = 270°) created by 
bifurcation of the equilibrium at circular orbit are referred to as 
Kozai equilibria, by analogy to the restricted problem (Kozai 1962; 
Lidov 1962). This change of stability of the central equilibrium in- 
duces large variation in eccentricity for an inner body initially on a 
nearly circular orbit, since its real motion (short periods included) 
will stay close to the separatrix of the reduced problem. The max- 
imal mutual inclination corresponding to the change of stability of 
the central equilibrium, called critical mutual inclination, has been 
calculated by iKozail j 19621 ): it drops from the well-known value 
39.23° to 32°, as the semi-major axes ratio increases from to 0.5. 
For the parameters of Figure [T] the critical mutual inclination is 
39.1°. 

Additional bifurcations of these equilibria occur for higher 
values of mutual inclination (see middle right panel of Figure [TJ. 
For more detail, we refer to t he complete study of these equilibria 
and their stability realized by IMigaszewski & Gozdziewskil J2009l) 
for the three-body problem. For increasing values of AMD, the 
equilibrium at the origin vanishes and the region of permitted mo- 
tion is divided into two islands. The dynamics is then governed 
by two families of equilibria: the equilibria related to the bifurca- 
tion of the Kozai equilibria and loc ated at the border of permit- 
ted motion (called solutions IVa by Mig aszewski & G ozdziewskil 
l2009i) ). and those related to the bifuracion of the central equilib- 
rium and appearing close to the e\ = Q axis (called solution Ilia 
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Figure 1. Level curves of constant Hamiltonian U9t in the representative plane (ei cos2gi,f2 cos Atti) for different values of AMD such that the mutual 
inclination at the origin is 20° (top left), 50° (top right), 80° (middle left), 110° (middle right), 180° (bottom left). In the bottom right panel, the origin 
ei = £2 = does not belong to the region of permitted motion. Other parameters are a = 0.05 and mi/m2 = 10"'*. 



bv lMigaszewski & Gozdziewskil ([20091)). Let us note that these last 
ones are unstable. 

In the Laplace plane reference frame, iLibert & HenrardI 

( 120081) have shown that, when the orbit is outside the Kozai- 
resonant area, the global ex trema of the eccentricities are reached 
when sinAw = (see also i Michtchenko et alj|2006r) , while their 
local extrema are reached when sin 2gi = 0. Example of such a be- 
havior is illustrated in Fig[T](top right panel) where the four points 
of intersection of a given orbit are symbolized by ' *' signs. As they 
are located in the four quadrants, the dynamics of this orbit is char- 



acterized by the circulation of both angles 2gi and Ao;. For a Kozai- 
resonant system conside red in the Laplace plane reference frame, 
iLibert & TsiganisI ( 120091) have shown that the eccentricities of both 
planets are not coupled, the eccentricity of the inner planet being 
extremal when sin 2gi = 0, and the one of the outer planet when 
sin Atij = 0. The ' + ' signs in Fig[T](top right panel) show a Kozai- 
resonant orbit: all the intersection points are located on the left part 
of the representation, indicating the libration of gi resonant angle. 
For the reasons given above, a particular interest of such a 
geometric view of the dynamics is to give an estimation of the 
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Figure 2. Level curves of constant Hamiltonian \\9\ in the representative plane (e\ cos2gi,e2 cos A^)) for different values of AMD such that the mutual 
inclination of the orbits with initial eccentricities ei = and £2 =0.1 is 20° (left), 50° (middle) and 80° (right). Other parameters are a = 0.05 and 
m[/m2 = 10""*. 



variation in eccentricity of each body. Let us note that, due to our 
choice of mass ratio fi, the eccentricity of the outer massive body 
is only weakly affected by its small inner companion. Indeed, the 
long-term variation in eccentricity is described by the Hamiltonian 
equation ( I20l l: 



V 



■«?a7C 



L,e, 



(23) 



Thus the variation of the outer eccentricity is of the order of mi, 
which is very small in this work, while the variation of the small 
body's eccentricity is quite important, as it is of the order of m^. 
As a result, the eccentricity of the perturbing body is nearly con- 
stant and it explains that all the intersection points of an orbit seem 
to have the same ordinate (in absolute value) in Fig. [T] (top right 
panel). On the contrary, the variation in eccentricity of the inner 
body can be very significant. For instance, the eccentricity of the 
orbit denoted by '*' in Fig [T] (top right panel) varies roughly from 
0.28 (positive abscissa) to 0.5 (negative abscissa). 

The representations of Figure [T] give information on the dy- 
namics of a system with nearly circular orbit of the perturber (e? 
smaller than 0.01). In the next section, dynamics with higher initial 
values of 62 will be considered. 

3.2 Elliptic orbit of the perturber 

In this section we assume an elliptic orbit of the perturber and study 
the dynamical evolution of the system by means of the representa- 
tive plane. While the dynamics of Figure[T]focus on nearly circular 
orbit of 1112, we keep increasing the AMD values to reach higher ec- 
centricities of this body. First let us consider £2 close to 0.1. Three 
different values of AMD are displayed in Figure |2] they are chosen 
such that the mutual inclination of the orbits with initial eccentric- 
ities ei = and £2 = 0.1 is 20° (left), 50° (middle) and 80° (right). 



Our first observation is the similarity to the dynamics of the bot- 
tom right panel of Figure[T] So the 3-D elliptic three-body problem 
is affected by two kinds of equilibria only: the "Kozai-bifurcated" 
equilibria at very high value of ei and the equilibria at circular orbit 
of the inner body. 

Although the dynamics is very similar for the three values of 
AMD displayed in Figure |2l the shifting on the y-axis is obvious 
and explains the different dynamics observed for a given system 
considered at various mutual inclination s. In order to explain ana- 
lytically the results of iFunk et al.l ( 1201 l i) (behavior of an Earth-like 
body initially on an inner circular orbit in a gas-giant system), let us 
consider the evolution of a two-planet system whose initial eccen- 
tricities are ei = 10"' and £2 = 0.1. The intersection points of the 
evolution of the system with the representative plane are denoted 
by '*' signs in Figure |2] 

For a small mutual inclination (i = 20°), gi and Ao) circulate 
and both variations in ei and 62 are so limited that the four expected 
intersection points seem gathered on two points only (see left panel 
of Figure |2j. For a mutual inclination of the orbits of 50° (mid- 
dle panel), the system is destabilized by the unstable equilibria: gi 
vacillates between libration and circulation and ej reaches values as 
high as 0.55. The same instability is present for the orbit of the third 
panel of Figure|2](i = 80°), where ei reaches a value close to 1. If 
we had extended the integration time, the intersection points would 
not be regular anymore, showing the chaotic evolution of highly in- 
clined systems due to the closeness to the "separatrix". These two 
different long-term evolutions (; = 20° and i = 50°) are illustrated 
in Figure [3] by means of a numerical integration of the octupole 
Hamiltonian equations J20I ). 

This change of dynamical behavior can be easily deduced 
from the shape of the Hamiltonian curves on the representative 
plane, as it can be seen on Figure |4] For small mutual inclina- 
tion (i.e. close to the borders of higher absolute values of ei), the 
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Figure 3. Long-term evolution of a system consisting of a small body initially on nearly circular orbit {e\ = 10" ) and an outer body whose initial eccentricity 
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systems whose inner orbit is circular correspond to extrema of the 
Hamihonian curve (left panel of Figure |4). On the other hand, for 
mutual inclinations higher than a value close to 40°, there exists 
another intersection point, belonging to the same curve of constant 
Hamiltonian, and of same eccentricity 62 (right panel of Figure |4j; 
the abscissa of these intersection points represent the secular varia- 
tion of the eccentricity of mi . For a higher eccentricity of the outer 
body, the dynamics is similar, as shown in Figure|5](e2 = 0.5). 

So we conclude that an inner small body on a quasi-circular 
orbit attracted by a giant companion on an elliptic orbit behaves 
secularly in a similar way as in the circular three-body problem: 
small periodic variation of its eccentricity when the mutual incli- 
nation of the orbital planes is small, on the contrary to the large 
chaotic variation observed for mutual inclinations higher than a 
critical value of ~ 40° . These ana l ytical results are consistent with 
the numerical study of lFunk et al.l ( 120 HI) . 

However, it is interesting to note that, even if the represen- 
tative planes of Figures [T] |2] and |5] precisely depict the dynamics 
around the central and Kozai families of equilibria, some addi- 
tional dynamical features can be "hidden". For instance, the 'O' 
symbols in Figure [T] (top right panel) identify an orbit character- 
ized by a libration of the angle hio around 180° (and simultaneous 
cir culation of gi). Th is kind of behavior is classified as mode 2 
bv lMichtchenko et al.l ( 1200 6) and its existence can not be deduced 
from the analysis of our representative plane. In the same way, no 
particular dynamics associated t o ~ 35 ° of mutual inclination, as 
the one reported by Fu nk et aj.! ( 120 111) and described in the next 
section, is visible on the representative planes of Figure [2] 



4 INTERESTING DYNAMICS AROUND 35° OF MUTUAL 
INCLINATION 

iFunk et alj j201lh have realized a numerical study of the long-term 



stability of inclined fictitious Earth-mass planets moving under the 
attraction of an eccentric giant planet. The small body is initially on 
a quasi-circular orbit. Although the massless companion exhibits 
periodically and limited variation in eccentricity for a mutual incli- 
nation smaller than ~ 40° (observation in agreement with our ana- 
lytic study of the previous section), their simulations have identified 
a dynamical region around 35° of mutual inclination consisting of 
long-time stable and particularly low eccentric orbits. 

This feature is well described by our octupole approxima- 
tion, as it is shown in Figure |6] where numerical integrations of the 
Hamiltonian equations l |20b are used to deduce the maximal eccen- 
tricity variation of the inner body on initially quasi-circular orbit 
(fii = 10"^), and for mutual inclination of the orbital planes up to 
50°. Left panel shows that, for small values of /, the variation of 
ei is negligeable when the orbit of the perturber is quasi-circular 
(due to the presence of the central equilibrium of Figure [T], while, 
for eccentric orbit of the perturber, the higher the value of ei the 
wider the secular variation of ei. This variation is even wider for 
high semi-major axes ratios, as shown in the right panel of Figure 
[6] For mutual inclinations higher than ~ 38°, the instability related 
to the Kozai bifurcations and described in the previous section pro- 
duces the important increase of the secular variations of ei (see the 
right sides of the graphs). The main new feature of Figure |6]is the 
sudden decrease of the maximal Si variation around a value close 
to 35° for all semi-major axes ratios and e ccentricities of th e gas 
giant, as observed in the numerical study of iFunk et al.l ( 1201 ll) . Let 
us note that this region around 35° is not present for a circular orbit 
of the perturber and is much in evidence for high semi-major axes 
ratios, so that the change of dynamics is due to the third-order terms 
of the octupole expansion. The aforementioned behavior is also il- 
lustrated, in Figure |7] by means of long-term evolutions of systems 
with mutual inclination of 34°, 35.55° and 38° (numerical integra- 
tion of the octupole Hamiltonian equations (I20t). In the following, 
this feature is analyzed in more detail. 
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To understand this particular behavior, we decide to realize an 
an alytical study of the frequ encies of the system, similar to the one 
of iLibert & HenrardI ( 1200 8h . Using a 12th-order expansion of the 
perturbative potential in powers of the eccentricities and the incli- 
nations, they have performed Lie transformations to introduce an 
action-angle formulation of the Hamiltonian and identify the an- 
alytical expressions of the four fundamental frequencies of the 3- 
D secular (non-resonant) three-body problem. This study has been 
realized in two reference frames: a general one and the Laplace 
plane reference frame. Our aim in the present section is to wonder 
whether the dynamics around 35° pointed out in Figures |6] and [7] 
corresponds to a commensurability between the fundamental fre- 
quencies. For the sake of completness, their analytical study is 
briefly described here. 

The Hamiltonian function expanded in powers of the eccen- 
tricities and the inclinations and averaged over the mean anomalies 
Mi reads: 



7C = - 



Gm\m2 



z 



D*-J1 J2 pljl l+2'l E-l;2l+2'2 c l*+;i l+2'l r. l*+J2l+2/4 



I'-JlJl-ilJe* 



with <D = [jipi - J2P2 -(k + ji)qi + {k + J2)qi\, Ef = y/IFjLi 
and 5, = s/lQj/Lj. The canonical variables in foimula ( 124) are the 
classical modified Delaunay's elements: 



/I, = mean longitude L; = m, yGmoo^ 

Pi = - longitude of the pericenter F, = L, 1 - -^1 - e,- 

g, = - longitude of the node 2, = L, .J 1 - e^ [1 - cos ;';] . 

(25) 
The indices (k, ii, / 6 4) are positive integers. The coefficients 
2j J1J2 (jgpgjj(] only on the ratio ai/oT of the semi-major axes. The 
secular Hamiltonian is a four degree of freedom problem. Let us 
note that it only depends on three angles, as 

* = ,/i(Pi -qi)- hiP2 - qi) - k(qi - q2). (26) 

As shown in lLibert & HenrardI (l2007h . the numerical convergence 
of the secular series (I24t is very good for a large set of parameters, 
even for moderate values of the eccentricities and the inclinations. 
The development is limited to order 12 in the eccentricities and 
the inclinations, which means that are kept in the Hamiltonian the 
terms for which the sum of the exponents of £1, £2, Si and 52 is 
lower or equal to 12. 

In order to obtain the analytic expressions of the four funda- 
mental frequencies, they have used a Lie transform perturbation 
scheme (Deprit 1969) to average the Hamiltonian ( 124) over the 
secular variab les p'. and q'. (i .e. the secular variables after a "reduc- 
ing rotation" ( iHenrardlll988h ). After a second Lie transform on the 
combination p\ + p'2 + q'l - iq'2, they get the following action-angle 
formu lation of the Hamiltonian ~ we refer to iLibert & HenrardI 
( I2OOS) for more detail-: 



■7C' 



z 



C, 



hJiJi 



y, 2/, p, 2/2 „-, 2/3 



(27) 



/l+/2+/3«12 



The associated Hamiltonian equations lead to the expression of the 
four frequencies: 

^_ (!-'') V9;/- p/2('l-l)p,2/2„-,2;3 



,-,2/3 



liJe3 

1^2 - ~lJ^2^^hCi,j2j,E[ £2 ' '^'1 

/,J£3 
^_ (!-/') V??/- r/2'lr/2'2o-,2(/3-l) 

/„/e3 

9l = 0, 



(28) 



Table 1. Long-term evolution of a system with r = 30° , obtained by decom- 
positions of frequencies on the data sets of the octupole approximation. Pe- 
riods are expressed in years. Initial pai'ameters of the system are f 1 = 10"*, 
£2 = 0.3, a = 0.05 and m\ /m2 = 10" 
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COS (i>, /J being the mass ratio m\l{mi + m2). The unit of frequency is the 
keplerian frequency 112 = JCmo/a^ of the mass OT2 multiplied by 
the mass ratio (mi -1- m^) I niQ. Let us note that, in the Laplace plane 
reference frame, the long-term evolution of the orbital elements can 
be described by only two frequencies and their linear combinations: 
/i = -P'l + q\ and /2 = -p'-, + q\ . 

In the following, we study the evolution of these frequencies 
with increasing values of the mutual inclination b etween the or- 
bital planes. By resorting to a frequency analysis (lLaskanll993h 
on the data sets obtained with the octupole approximation. Table [T] 
identifies the main combinations of the fundamental frequencies 
common to the evolutions of the orbital elements for i = 30°. The 
frequencies are listed by decreasing amplitude of the trigonometric 
term and noted ci (highest amplitude) to C5. Bold type ci corre- 
sponds to the precession rate of an angular variable in circulation. 
Last columns display the identifications of the different combina- 
tions in terms of the fundamental frequencies (p[ , p'2 and q[ ) and 
the two frequencies /i and /2 respectively. 

As can be observed in Table \T\ the two frequencies, /i = 
0.01184 and /2 = 0.00837 (values calculated from equation ^ 
of the analytical I2th-order expansion, in their unit of frequency), 
correspond to the precession rates of the arguments of the pericen- 
ter (1)1 and ii)2 respectively. Last column shows that all the frequen- 
cies of the orbital evolutions are linear combinations of these two 
frequencies. In particular, the main frequency of the eccentricities 
is the precession rate of Aw and corresponds to /i - f2- Let us note 
that the analytical frequencies given by equation ( 128) are very close 
to the ones identified by the frequency analysis on the octupole ap- 
proximation: fiaci = 0.01192 and f2oct = 0.00839. This small shift 
in the periods (less than lO' years) is due to the limitations of both 
approximations with respect to the semi-major axes ratio or the ec- 
centricities and inclinations. 

To analyze the dynamics around a mutual inclination of 35° 
in more detail, we examine the evolution of the two frequencies for 
increasing mutual inclination values. Figure |8]shows the evolution 
of the periods associated to /i and /2 for mutual inclination up to 
38°. The two curves intersect when i =~ 36°, namely 35.55° in the 
octupole formulation and 35.9° in the eccentricities and inclinations 
development. 

As a result, the coupling between the frequencies of the orbital 
elements is different for this particular value of the mutual incli- 
nation, as shown in Table [2] (frequency analysis on the octupole 
approximation). Indeed, all the frequencies are combinations of 
f = fi = /2 and a very small frequency g. The change of dy- 
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Table 2. Same as Table[T]for i = 35.55° 
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Figure 9. Long-term evolution of a system with mutual inclination of 
35.55° in a frame where the reference plane is the plane of the giant planet. 
Note that the resonant angle is (ii\ -(1\. 



namics induced by the commensurability f\ = fi, is obvious when 
looking at the Aco's evolution in Figure |7] the angle is in libration 
for the very low eccentric orbit at / = 35.55°. As it corresponds to 
a behavior modification of an angle, this particular dynamics can 
thus be regarded as a secular resonance. 

All this study is realized in the Laplace plane reference frame. 
If we consider the orbital plane of the giant planet as reference 
plane ((2 = 0), the evolution of the eccentricities is similar and the 
resonant angle becomes aii - fij, as illustrated in Figure |9] 



5 CONCLUSION 

In the present work, we focused on the study of the 3-D elliptic 
three-body problem of a small mass under the attraction of an outer 
giant body. The influence of an eccentric orbit of the perturber on 
the dynamics of a small inclined inner body has, to our knowledge, 
not yet been investigated in the literature. Particular attention has 
been given to a region around 3 5° of mutual inclination detected 
numerically bv lFunk et al.l ( 120111) . 

Our analytical study relies on the octupole expansion, which 
is a compact formulation of the Hamiltonian suitable for hierarchi- 
cal planetary systems. Short-period averaging and node reduction 
(by adoption of the Laplace plane reference frame) enable us to re- 
duce the problem to two degrees of freedom. The four-dimensional 
dynamics is analyzed through representative planes which identify 
the main equilibria of the problem. It has been shown that an inner 
body on quasi-circular orbit behaves secularly in a similar way as in 
the circular three-body problem: its eccentricity variations are very 



limited for mutual inclination between the orbital planes smaller 
than ~ 40°, while they become large and chaotic for higher mutual 
inclination. 

As shown by iFunk et alj ( 1201 11) . there exists a dynamical re- 
gion around 35° of mutual inclination consisting of long-time sta- 
ble and particularly low eccentric orbits of the small body. Using 
a 12th-order Hamiltonian expansion in eccentricities and inclina- 
tions, in parti cular its actio n-angle formulation obtained by Lie 
transforms in iLibert & Hen rard (2008), we have shown that this 
region corresponds to a commensurability of the two frequencies 
that are the precession rates of the arguments of the pericenter oi, 
and (i)2- It explains the change of dynamics of the angle Aw which 
starts to evolve in libration. This particular dynamics can thus be 
regarded as a secular resonance. The same analysis can be realized 
with the orbital plane of the giant planet as the reference plane (i.e. 
no adoption of the Laplace plane) to identify a>i -Q.^ as the resonant 
angle of this reference frame. 

This study also applies to binary star systems, where a planet 
is revolving around one of the two stars (inner problem), since the 
mass ratio between the bodies is of the same order (/j ~ 10"**). 

The region around 35° could belong to the habitable zone of 
the system and be of particular interest for the research of life in 
extrasolar systems, as it consists of stable orbits with limited ec- 
centricity variation of the planet, which means a constant distance 
between the planet and the host star. 
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